NASA Technical Memorandum 100218 

ICOMP-87-6 

AIAA-88-0594 


Navier-Stokes Cascade Analysis With 
a Stiff k-e Turbulence Solver 

N 8 8- 10778 

Unclas 
0106497 

Jong-Shang Liu 

Institute for Computational Mechanics in Propulsion 
Lewis Research Center 
Cleveland, Ohio 

Peter M. Sockol 
Lewis Research Center 
Cleveland, Ohio 

and 

Joseph M. Prahl 

Case Western Reserve University 

Cleveland, Ohio 


(NAS A-TN- 100 21 8) NAVIEB-STOKES CASCADE 
ANALYSIS WITH A STIFF KAPPA -E PSILON 
TURBULENCE SOLVER (NASA) 15 p Avail: NTIS 
HC A03/HF A01 CSCL OtA 

G3/02 


Prepared for the 

26th Aerospace Sciences Meeting 
sponsored by the American Institute of 
Reno, Nevada, January 11-14, 1988 




Aeronautics and Astronautics 



NAVIER-STOKES CASCADE ANALYSIS WITH 
A STIFF k-e TURBULENCE SOLVER 


Jong-Shang-Liu* 

Institute for Computational Mechanics in Propulsion 
NASA Lewis Research Center, Cleveland, Ohio 

Peter M. Sockol 

Computational Fluid Dynamics Branch 
NASA Lewis Research Center, Cleveland, Ohio 

and 


Joseph M. Prahl 

Department of Mechanical & Aerospace Engineering 
Case Western Reserve University, Cleveland, Ohio 


Abstract 

The two-dimensional, compressible, thin-layer Navier- 
Stokes equations with the Baldwin-Lomax turbulence 
model and the k-€ model axe solved numerically to simu- 
late the flow through a cascade. The governing equations 
are solved for the entire flow domain, without the bound- 
ary layer assumptions. The stiffness of the k-£ equations is 
discussed. A semi-implicit, Runge-Kutta, time- marching 
scheme is developed to solve the k-e equations. The im- 
pact of the k-e solver on the explicit Runge-Kutta Navier- 
Stokes solver is discussed. Numerical solutions are pre- 
sented for two-dimensional turbulent flow over a flat plate 
and a double circular arc cascade and compared with ex- 
perimental data. 

Introduction 

In order to reduce the number of trial designs of tur- 
bomachinery, improved understanding of the physical phe- 
nomena occurring within the turbomachine is required. 
A wide variety of flow characteristics that exist in tur- 
bomachinery blade passages is described in detail by Mc- 
Nally and Sockol 1 : large pressure gradients, shock waves, 
shock-boundary layer interactions, and cross flows gen- 
erated by vorticity field redistribution in turning flows. 
With the advent of large digital computers and the devel- 
opment of efficient numerical algorithms, numerical simu- 
lation models have become a realistic alternative to labo- 
ratory testing programs. The flow field predicted by nu- 
merical calculations can be validated with available data 
from experiments. If acceptable agreement is obtained, 
a greater level of confidence in the predictive capability 
of the numerical model will enhance the usefulness of the 
code to the designers. 

The flow within turbomachinery may be laminar, 
transitional or turbulent, attached or separated, and 
wakes exist downstream of the blade rows. To obtain 
useful numerical predictions for such complex flows, the 
Navier-Stokes equations are required. Although the 
Navier-Sfcokes equations are valid for turbulent flows, the 
fluid variables change rapidly in space and time and 
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can not be resolved by numerical calculations. After the 
compressible form of the Navier-Stokes equations are mass 
averaged, a turbulence model is required to obtain a closed 
system of equations. 

The purpose of this work is to apply two turbulence 
models to a cascade and compare the results with avail- 
able data. The two dimensional flow is studied by nu- 
merically solving the mass-averaged, compressible, thin- 
layer Navier-Stokes equations in conjunction with an al- 
gebraic eddy viscosity and two-equation model. The al- 
gebraic eddy viscosity model was developed by Prandtl 
in 1925, and modified by Cebeci et al. 2 as a two layer 
model. Baldwin and Lomax 3 modified it to avoid finding 
the edge of the boundary layer. The two-equation model, 
developed by Jones and Launder 4 , employs less empiri- 
cism than the algebraic eddy viscosity model: the veloc- 
ity and length scale are calculated. It contains two addi- 
tional partial differential equations, the turbulence kinetic 
energy (k) and dissipation rate (e) equations. The k and e 
equations used here are in a compressible form without the 
boundary layer type assumptions. A two dimensional vis- 
cous code developed by Chima 5,6 using a finite-difference, 
explicit, two-stage, time-marching scheme is used to pre- 
dict the mean flow. Since the source terms of the k and 
£ equations are stiff, a semi-implicit, Runge-Kutta, time- 
marching scheme is developed to calculate the k and € 
equations. 

The turbulent flow over a flat plate and the flow 
through a two dimensional double circular arc cascade are 
calculated. The comparison between the calculated and 
measured flow fields are presented and discussed for both 
the Baldwin-Lomax method (BLM) and the Jones and 
Launder k-6 method (KEM). 

Formulation 

Governing Equations 

The equations governing the fluid motions are the 
compressible form of the Navier-Stokes equations. For 
turbulent flow, the equations are mass averaged. The 
system of equations is solved numerically on a body- 
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conforming grid system, so the equations are transformed 
from the Cartesian system (x, y) to an arbitrary body- 
fitted coordinate system (f, 77 ). The thin-layer approxi- 
mation is used in the transformation. It neglects the vis- 
cous diffusion terms parallel to the body surface. 

The equations are nondimensionalized by reference 
quantities, Chima 5 , as follows, 

i=t/(L*/c*), x = x/L*, y = y/L\ u = u/c*, 

V = vjc, p = p/(pV 2 ), f = T/To, e = e/( P V 2 ), 

H = H/c t2 , p = p/p\ fi = p/n*, 

D P*c*L* D 

— — , Re= — , Pr =~~ * 

'v 4- 1 . u, A* 

The nondimensionalized, conservative form, two- 
dimensional, time-dependent, mass averaged, thin-layer 
Navier-Stokes equations are written as follows, (neglecting 
the overbar) 


Two-Equation Model (KEM) 

In this model, first introduced by Jones and Launder 4 , 
the eddy viscosity is obtained from the turbulence kinetic 
energy and dissipation rate. The turbulence quantities, k 
and c are determined from the transport equations. 

The k and e are nondimensionalized by k* and £*, 
where k* = ii*c m /{p*L*), e* = fi*c* 2 /(p* L* 2 ). The 
nondimensionalized, thin-layer form of the k-e equations 
are written in a conservative form similar to equation (l). 

*-( 2 ) *-(#) *-(£) 

r _ ((P + Pt/°k)™ sk„\ r — ( S ~P e \ 
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where 

mi = {4/3t} 2 -f rj 2 )(fi + pt) 

m 2 = l/3ri x riy{tJi + fit) 
m 3 = (t7 x + 4/3rjl)(fi + iit) 
m 4 = {nl + ril){n/Pr + Pt/Pr t ) 


Turbulence Model 


Baldwin-Lomax Model (BLM) 

This model is the two layer model developed by 
Baldwin and Lomax 3 , and has been extensively used in the 
literature for calculating attached and mildly separated 
turbulent flows. 

The Prandtl-Van Driest formulation is used in the 
inner region to obtain pL t . A modified Clauser formulation 
is used in the outer region. A simple transition criterion 
is formulated in this model by setting the eddy viscosity 
equal to zero everywhere in a profile for which the 
maximum tentatively computed value of is less than a 
specified value. The detailed formulation is omitted here 
and can be seen in ref. 3. 


m 5 -n\ + rjl, fi t = C^pP/e, 
ni =4/3 £ + i 2 y , n 2 = 4/3 nl + 

n 3 “£x + 4 /3£ 2 , n 4 - 7? x + 4/3?7^, 

n$ =8/3£ x 77 x + 2^y7/y, tiq = 8/3£y7/y -f* 2 £ X T 7 X , 

n 7 =2/ 3£ x fy, tiq ~ 2/377 x 77y, 

ng = 2 £ x ? 7 y — 4 / 3 ^y 7 / x , 7110 = 2 £yf/ x — 4 / 3 £ x 77 y, 

5 — + 712 ^^ + 713 ^ + n 4 vl + nsu^Urj + 

+ n 7 U£V£ + neu v vrj + ngu^v^ + tiiqu^v^) 

- 2/3pk(u^ x + U v nx + V$£y + Vr,n y ) 


According to Launder and Spalding 8 , the constants 
appearing in (3) are <7*=1, a £ =1.3, Ci=1.44, C2=l-92, 
=0.09, and no effort is made to change these numbers 
in this study. 

The modeled k and e equations are effective only in 
the fully turbulent region. 9 Treatment of the near wall 
region is described in the boundary condition section. 

Numerical Scheme 


A two dimensional viscous code developed by Chima 5 ’ 6 
is used to predict the mean flow. The code employs an 
explicit, finite-difference, two-stage, Runge-Kutta algo- 
rithm. 


qi = F< n] - a, AtR{F< n) ) 
q 2 = F 0 (n) - A tR( qi ) 


= 1.2 


R =\\(Fii +1 - hly) + (F 2 i +l - F,’- 1 ) 
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The second and fourth order artificial dissipation terms 
are added to stabilize the scheme. 

For the k and e equations, it is unsuitable to use an 
explicit method, since the source terms of these equations 
are stiff. The integration of the stiff equations using an 
explicit method often imposes severe At limitations and 
has stability problems. A semi-implicit treatment of the 
source term of the k and e equations is used. A two-stage, 
semi- imp lie it, Runge-Kutta scheme developed by Liu 7 is 
used to calculate the k and e equations. 

</i = A t[I - A<a J 4(F 0 (n) )l- 1 i2(F 0 (n) ) 

9 2 = A t[I - A<aA( J F 0 (,,) )]- 1 J R(F 0 (n) + b qi ) (5) 
^ +1) = H n) + a l9l + <22*72 

R =|[(Af +1 - ALi) + (£4 +1 - hr 1 ) 

A =dF 3 /dF 0 

The coefficients a, b, ai, and 0:2 are obtained by matching 
(5) with the Taylor series expansion of Fq H+1 ^ to second 
order in time and by letting the characteristic root go to 
zero as t goes to infinity. 

a = b = <*! = 1.7071, a 2 = -1.7071 

The convection terms are evaluated explicitly and the 
source term semi-implicitly. The two by two Jacobian 
matrix from the source term is evaluated once for every 
two stage calculation. A local time step is used to accel- 
erate convergence to a steady state. 

The impact of the k-e equation on the main code is 
introduced in the following way. 

(i) The main flow field is calculated with the BLM 
method for 1000 time steps. 

(ii) The flow field variables are fixed and the k and e 
equations are calculated for 1000 time steps. 

(iii) The k and e values are fixed and the mean flow field 
is calculated for 10 time steps. 

(iv) The flow field is fixed and k and e are calculated for 
10 time steps. 

(v) Step (iii) and (iv) are repeated until convergence to 
steady state is achieved. 

Stability Analysis 

The numerical stability is analyzed in the inviscid 
parts of the governing equations. By linearizing the in- 
viscid parts of the conservative form of the equations and 
applying the Von Neumann stability criterion, the stabil- 
ity limitation on the time step for the mean flow and k-e 
equation can be found. The At limitation for the mean 
flow is 

A t < 1/[U + V + cy/(£ t + r/x) 2 + (£y + r/y) 2 ]. 


The At limitation for the k-e equation can not be shown 
in explicit form, see ref. 7 for the detail. 

Boundary Conditions 

For subsonic inflow, the three incoming characteris- 
tic speeds are replaced by three boundary conditions: 
total pressure, total temperature, and flow angle. Small 
values of k and e are specified at the inlet. The upstream- 
running Riemann invariant is extrapolated from the inte- 
rior. For subsonic outflow, the static pressure is specified 
at the downstream boundary and other properties are ex- 
traploated from the interior. No-slip conditions and the 
wall temperature are specified on the wall. The pressure 
on the wall is calculated through the normal momentum 
equation. 

Since the k and e equations are only effective in a 
fully turbulent region, they are not applied directly to the 
wall. Many modifications are required to calculate k and 
6 near the wall. A review of the low Reynolds number 
two-equation model given by Patel, Rodi and Scheuerer 9 
bases the modifications largely on numerical experiments 
and comparisons with different flows, and it is not clear 
which models can be used with confidence. Near the wall, 
the values of k and e change very rapidly, and many grid 
points are needed to predict k and e there. In this study, 
the BLM method is used to predict in, and from the 
assumption of local equilibrium, values of k and e near 
the wall are found, (y + < 40). This method does not 
need to modify the equations, is easy to calculate, and 
does not appear to require very fine grids near the wall. 

Results and Discussion 

Numerical results are presented for turbulent flow 
over a flat plate and for flow over a two-dimensional double 
circular arc cascade. 

Flat Plate Problem 

The subsonic flow over a zero thickness flat plate 
cascade with zero stagger angle is calculated comparing 
the KEM and BLM turbulence models. The flat plate of 
chord C and zero thickness has a pitch of 0.2C. Because 
of symmetry, the computational domain is one half of the 
flow passage. The inlet boundary is 0.1C upstream of the 
leading edge of the plate. The flow with an inlet Mach 
number of 0.5 and zero angle of attack is calculated for 
this case. The minimum grid space Ax is 0.0025C in the 
near leading edge region and smoothly stretched up and 
downstream. The minimum grid space Ay is 10“ 6 C and 
the grid is smoothly stretched away from the wall. 

In Fig. 1 the computed local skin friction coefficient 
versus Reynolds number based on the distance from the 
leading edge is compared with the Karman-Schoenherr 
skin friction correlation 10 , Wieghardt 11 data and the 
correlation from Schlichting 12 . The numerical results of 
the two models compare well with the data, with the 
BLM predicting a higher Cf than that predicted by KEM 
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for Re x greater than I0 6 . The predicted behavior in the 
range 10 5 < Re x < 10 6 is due either to the effect of 
the singular point at the leading edge or to the flow not 
being fully turbulent in that region. Figs. 2 and 3 show 
the velocity distribution. In Fig. 2, the mean velocity, 
nondimensionalized by the free stream velocity, is plotted 
versus the distance from the wall, nondimensionalized 
by the displacement thickness. KEM and BLM predict 
the same results in excellent agreement with the data of 
Klebanoff and Diehl 13 . In the near wall region, u + is 
shown versus y + in Fig. 3. The numerical results agree in 
the sublayer region. In the log-wall region, the predicted 
Karman constant is about 0.4, however the predicted u + 
is smaller than data of Schlichting 12 , Klebanoff 13 and 
Hinze 14 . This is most likely due to a lack of enough grid 
points in the region 10 < y + < 50. 

The distribution of the turbulent shear stress at Res 
= 80000 is compared with the data of Klebanoff 15 in Fig. 
4 and that of Schubauer 16 in the near wall region in Fig. 
5. Near the wall, the computed turbulent shear stress is 
nearly constant and agrees favorably with the data. The 
discrepancies increase as y/S increases. This is probably 
attributable to grid size effect. The BLM predicts slightly 
higher turbulent shear stresses than that predicted by 
KEM. The KEM and BLM results are virtually identical 
in the near wall region and agree within experimental error 
with the data (Fig. 5). 

The turbulent kinetic energy and dissipation rate 
profile from the KEM calculations are shown in Figs. 6 
to 8. The computed k distribution is compared to the 
data of Klebanoff 15 in Fig. 6. The computed k agrees well 
with the data except near the wall. The k distribution 
in inner variables is compared to the results of Patel 
et al. 9 in Fig. 7. k + is underpredicted in the inner 
region and no peak of k + is predicted because the KEM 
is only effective in the fully turbulent region. In the 
sublayer and buffer region (y + < 40), the assumption 
that turbulence production is equal to dissipation rate 
is used. The dissipation rate is over-predicted in these 
regions (Fig. 8). Some modified versions of the two- 
equation model can predict the existance of the peak of 
k + . g However, most of the modifications are based on 
numerical experiments, without physical basis and applied 
to specific flows. Since the basic two-equation model can 
predict the mean flow well except very near the wall, the 
local equlibrium assumption and the zero-equation model 
are used in the near-wall region where KEM is ineffective 
or requires further questionable modifications. 

Double Circular Arc Cascade 

Numerical calculations for flow over a two dimen- 
sional double circular arc cascade are compared with data 
obtained on a blade section designed by Sanger at NASA 
Lewis Research Center and tested at Pennsylvania State 
University by Deutsch and Zierke 17 ’ 18 . 

The blade section is a double circular arc with a 
camber of 65 degrees and a chord length of 228.6 mm, 


machined round to a 9.14 ym radius at the leading 
and trailing edges. The equations for the blade surface 
geometry, and camber line are given in 17 . In the Penn, 
state experiment, a five blade cascade was tested. The 
solidity (chord length/pitch length) of the cascade was 
2.14 and the stagger angle was 20.5 degrees in the test. 
The measurements were made in air at a chord Reynolds 
number based on inlet conditions of 500,000 and a positive 
incidence angle of 5 degrees. The inlet velocity, 33.11 
m/sec (about Mach 0.1), was measured by a five-hole 
probe approximately 38 mm upstream of the blade leading 
edge. A one-component laser Doppler velocimeter (LDV) 
was used to measure the flow field in the experiment, 
and the five-hole probe was used to make inlet and outlet 
measurements. The two dimensionality in the experiment 
was controlled by adjusting the lower and upper false 
blade positions, the tailboard position, and the magnitude 
of top, side and upper-channel suction. The periodicity of 
the cascade was controlled by side wall suction. 17 

The computational domain covers one blade passage 
of the cascade. The upstream boundary is located about 
sixty percent chord upstream of the leading edge and 
the downstream boundary is located about sixty percent 
chord downstream of the trailing edge. The interior 
body-fitted grids of this domain are generated using 
the GRAPE code developed by Sorenson 19 . This grid 
generation code generates a two dimensional grid about 
arbitrary boundaries by iteratively solving two coupled 
Poisson equations. Two C-type grids with (257 x 65) and 
(129 x 33) points are used in the numerical calculations. 

Fig. 9 shows a comparison of the calculated and mea- 
sured blade static pressure distributions. The pressure 
surface has a large favorable pressure gradient near the 
leading edge. As the flow proceeds, the pressure gradi- 
ent is mildly adverse and becomes favorable again further 
downstream. The suction surface has a large unfavorable 
pressure gradient near the leading edge and a flat pressure 
profile near the trailing edge as seen from the data. The 
pressure distribution on the pressure surface predicted by 
the coarse grid calculation is identical with the fine grid 
calculation. The predicted pressure distribution on the 
suction surface by the fine and coarse grids are different 
before 40 percent cord and after 70 percent chord because 
the resolution of the coarse grid is not fine enough to pick 
up the rapidly changing flow properties on the suction sur- 
face. The results predicted by a laminar flow calculation 
also shown in Fig. 9 do not predict the turbulent flow well, 
especially on the suction surface. The KEM and BLM re- 
sults show similar results for pressure distribution. Both 
predictions are in good agreement with data over the for- 
ward portion of the blade. The real flow tends to separate 
and becomes unsteady after sixty percent chord on the 
suction surface, and the discrepancies increase especially 
on the suction surface. The KEM predicts slightly better 
agreement with data than BLM on the aft portion of the 
blade. 

Calculated pressure and Mach contour are given in 
Fig. 10 and 11. No measured contours are available. 
Large pressure gradients occur near the leading edge on 


4 



both the pressure and suction surface. Small pressure gra- 
dients exist on the pressure surface aft of the leading edge 
while still large pressure gradients persist on the suction 
surface. In the wake, there are small pressure gradients 
indicated on the pressure contours. The stagnation point 
is on the pressure surface because of the high incidence an- 
gle. The high velocity generated near the leading edge of 
the suction surface is a result of flow acceleration from the 
stagnation point around the blade leading edge. When 
the flow proceeds further on the suction surface, and is 
subjected to a large adverse pressure gradient, the flow 
becomes unstable and transition to turbulence happens 
rapidly. 17 The growth of the boundary layer on the suc- 
tion surface and the wake development downstream of the 
trailing edge are evident in the Mach number contours 
(Fig. 11). 

The velocity distributions were measured at eleven 
chordwise locations on both the pressure and suction sur- 
faces, and the data have ninty-five percent confidence 
bands, (except near the leading edge thin boundary layer 
and near the trailing edge unsteady flow). 17 Since the 
computational grids are non-orthogonal and not coinci- 
dent with the measurement locations, two-dimensional lin- 
ear interpolation 20 is used to obtain the computed veloc- 
ities at the measurement locations to compare with data. 
Comparisons of the calculated and measured velocity dis- 
tributions at six chordwise locations on the suction surface 
are shown in Fig. 12. y is the normal distance from the 
wall and u is the velocity component parallel to the wall. 
The flow over all of the suction surface is turbulent be- 
cause the large adverse pressure gradient near the leading 
edge triggers early transition. The calculations agree with 
the data well near the leading edge, with the BLM show- 
ing marginally better agreement with the data than the 
KEM. Both models predict less flat velocity profiles than 
the data, possibly because of deficiencies in the turbulence 
models which can not handle high curvature and pressure 
gradient effects well or because insufficient resolution of 
the boundary layer in the numerical calculations under- 
predicts the turbulence near the wall. The KEM results 
have better agreement than BLM after 50 percent chord. 
After 60 percent chord on the suction surface, the flow 
tends to separate and become unsteady. The KEM pre- 
dicts separation after 90 percent chord and no separation 
is predicted by the BLM calculation. 

Fig. 13 compares the calculated and measured veloc- 
ity distributions at six chordwise locations on the pres- 
sure surface. At the leading edge on the pressure surface, 
because of the large incidence angle, the favorable pres- 
sure gradient maintains a laminar flow. The data shows a 
mildly adverse pressure gradient following the leading edge 
high favorable pressure gradient. Further downstream, 
the pressure gradient again becomes favorable and a fully 
turbulent flow does not develop on the pressure surface. 
The results from both the KEM and BLM calculations are 
similar on the pressure surface. The flow on the pressure 
surface is well behaved with no separation since there are 
no strong adverse pressure gradients to trigger separation, 
transition, or reattachment. 


Measurements of the wake were made at three lo- 
cations, two near wake and one far wake locations. The 
near wake data were taken at 105.4 and 109.6 percent 
chord positions. The percent chord given here is actually 
the percent of total pressure surface blade arc. The origin 
of y is along the extended arc from the pressure surface 
and y is positive on the pressure side of the curve. The 
far wake data were taken at 152.6 percent chord position. 
The origin of y is located along the tangent curve of the 
camber line at the trailing edge. The LDV technique was 
used to measure the flow field in the near wake and the 
five-hole probes were used to measure the flow in the far 
wake. Fig. 14 shows the velocity distribution in the near 
wake. The flow in the near wake is asymmetric because 
the suction surface and pressure surface boundary layers 
that merge at the trailing edge are different in shape and 
thickness. The KEM results are closer to the data than 
the BLM results, but there remains substantial disagree- 
ment between prediction and measurements. The veloc- 
ity distribution data and predictions in the far wake shown 
in Fig. 15 are approaching symmetric profiles. The center 
line of the wake is moved to the suction side because of 
the smaller velocity near the trailing edge on the suction 
surface than on the pressure surface, and because of the 
pressure difference between these two surfaces. The com- 
puted center line of the far wake is closer to the pressure 
side than that measured because the computed trailing 
edge velocity on the suction surface is higher than that 
measured. 

Conclusion 

A numerical solution of the two dimensonal, com- 
pressible, turbulent cascade flow field problem comparing 
the Baldwin-Lomax (BLM) and the k-c (KEM) turbulence 
models is solved using a two-stage, explicit, Runge-Kutta, 
time-marching code for the mean flow. The stiffness of 
the k and e equations necessitates a semi-implicit, Runge- 
Kutta, time marching scheme to solve the thin-layer type 
k-£ equations. No special effort is made in modifying the 
model constants to improve agreement between predic- 
tions and measurements in this study. The local equilib- 
rium assumption and the BLM is used to predict the val- 
ues of k and c near the wall in the KEM calculation. No 
wall function or low Reynolds number k-e modification 
is used. Numerical solutions are obtained for two dimen- 
sional, compressible, turbulent flow over a flat plate and 
a double circular arc cascade and compared with data. 

Flow over a flat plate with an inlet Mach number 
of 0.5 is calculated comparing the KEM and BLM turbu- 
lence models. Both the KEM and BLM results are similar 
and compare favorably with experimental data. For this 
simple geometry with no pressure gradients, flow separa- 
tion or surface curvature, the BLM is sufficient to predict 
the mean flow field. The use of a higher-order turbulence 
model does not have advantage over the BLM. 

For the double circular arc cascade flow at an inlet 
Mach number of 0.1, both the KEM and BLM predict 
nearly identical surface pressure distributions. The KEM 


5 



predicts velocity profiles that agree better with data than 
those predicted by BLM near the trailing edge separa- 
tion region on the suction surface. Although the KEM 
results are closer to the data than the BLM results, there 
remains substantial disagreement between prediction and 
measurements. The discrepancies between the computa- 
tions and the experimental data arise from both inadequa- 
cies in the turbulence model and numerical effects. The 
grid resolution in the surface normal direction and in the 
streamwise direction have an influence on the numerical 
accuracy. An adapative grid may be needed in the cal- 
culation where the flow properties change rapidly. The 
thin-layer approximation in the calculation might be un- 
suitable since the flow properties change rapidly along the 
streamline, especially near the leading and trailing edges. 

For flow over complex geometries with pressure gra- 
dients or separation, the BLM is not adequate and a 
higher-order turbulence model is needed. The application 
of KEM to complex flow geometries needs improvements 
to handle the low-Reynolds number flow near the wall, 
flow transition, and separation. It also requires improve- 
ments in computational efficiency to make the KEM cal- 
culation competitive with the BLM calculation. 
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Fig. 1 Comparisons of the computed and measured 
local skin friction coefficient over a flat plate 



Fig. 2 Comparisons of the computed and measured 
velocity distribution 
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o Re, = 2000 closed Symhoh KEM 

A Re e = 6000 Open Symbol: BLM 

d Re 0 = 10000 




Fig. 3 


u + = 5.75log(y + ) + 4.2(/c = 0.4) [13] 
u + = 2.44/n(y + ) + 4.9(ic = 0.41) [14] 
u + = 5.75 log(y + ) + 5.5(/c = 0.4) [8] (Pipe flow) 



Comparisons of the computed and measured 
velocity distribution in inner variables 


— KEM Prediction at Re s « 80,000 



Fig. 4 Comparisons of the computed and measured 
turbulent shear stress distribution 
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Fig. 5 Comparisons of the computed and measured 
turbulent shear stress distribution in inner 
variables 
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Fig. 7 The turbulent kinetic energy distribution in 
inner variables 



Fig. 8 The dissipation rate distribution in inner 
variables 
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Fig. 9 Comparisons of the computed and measured 
blade static pressure distribution 


Pressure increments ( Ap/po ) = 0.001 



Fig. 10 Calculated pressure contours 


Mach number increments (AM) =0.005 



Fig. 11 Calculated Mach number contours 
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Fig. 12 Comparisons of the computed and measured 
velocity distribution on the suction surface 


— KEM Prediction 
— BLM Prediction 
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Fig. 13 Comparisons of the computed and measured 
velocity distribution on the pressure surface 
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— KEM Prediction 

— BLM Prediction 

o Exp Data [17], [18] 



Fig. 14 Comparisons of the computed and measured 
velocity distribution in the near wake 



Fig. 15 Comparisons of the computed and measured 
velocity distribution in the far wake 
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